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ABSTRACT 

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains 
the physics of cosmic rays (CRs) in these dense starbursts. We argue that the 7-rays are predominantly hadronic 
in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR 
protons in both galaxies: we show that a fraction F ca i « 0.2-0.4 of the energy injected in high energy primary 
CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing 7-rays, 
neutrinos, and secondary electrons and positrons. We discuss the factor ~ 2 uncertainties in this estimate, 
including supernova rate and leptonic contributions to the GeV- TeV emission. We argue that 7-ray data on 
ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits 
on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts 
suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. 
We briefly reconsider previous predictions in light of the 7-ray detections, including the starburst contribution 
to the 7-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest 
star-forming galaxies on the sky and present updated predictions for their 7-ray and neutrino fluxes. 
Subject headings: galaxies: individual (M82, NGC 253), starburst - cosmic rays - gamma rays: theory, obser- 
vations - radio continuum: galaxies 



1. INTRODUCTION 

M82 and NGC 253 are nearby (D « 2.5-4.0 Mpc), proto- 
typical starburst galaxies, each having an intense star-forming 
region of about 200 pc radius in the center of a more quies- 
cent disk galaxy. The starbursts are expected to have high 
supernova (SN) rates of about 0.03-0.3 yr . SN rem- 
nants are believed to accelerate primary cosmic ray (CRs) 
protons and electrons. The high SN rates in starbursts im- 
ply high CR emissivities. The presence of CR electrons and 
positrons in these starbursts is inferred from the nonthermal 
synch rotron radio emission they produce (e.g., iKlein et all 
119881: IVolket al.lll989h . However, most of the CR energy is 
believed to be in the form of CR protons. 

When high energy CR protons collide with interstellar 
medium (ISM) nucleons, they create pions, which decay 
into secondary electrons and positrons, 7-rays, and neu- 
trinos. With their dense ISMs ((«) « 100-500 cm" 3 ) 
and possible high CR en ergy densit i es (as evinced by 
the b r ight radio emissi on jVolk et al.1 [19891: lAkvuz et alJ 
[19911 iPersic & Rephaehl |2010|), M82 a nd NGC 253 are 
predicted to be bright 7-rav sources (e.g.. | Akvuz et al!l | 199 l b 
Sreekumar et1dll994IVolk et alj[~1 996: Pag lione et all 199a 
Romero & Torres 2003; Domingo-Santama rfa & Torres! 
20051: iThompson. Quataert. & Waxmanl 120071 [ TQW1 ; 

120091: 
As pro- 



Persic etal 



Rephaeli et al. 



Ide Cea del Pozo et alJ 
lLackietal.1 120101 rLTOl). 



totypical starbursts, if M82 and NGC 253 are seen in 7-rays, 
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starbursts in gener al may be sources of 7-rays (lPohllll994t 
iTorres et all 12004a) . with important implica tions for the dif- 
fuse 7-ray and neutrino backg rounds (e.g. JPavlidou & Fields! 
12001 lLoeb & Waxiiianll200a rLW061: TQW). However, the 
7-ray luminosity of starbursts depends not only on the injec- 
tion rate of CRs, but also on the efficiency of converting CR 
proton energy into pionic 7-rays, neutrinos, and secondary 
electrons and positrons. In turn, this efficiency depends on 
the ratio of the timescale for pion production to the escape 
timescale. The hypothesis that CR protons in starbursts 
lose all of their energy to pionic c ollisio ns before escaping 
is called "proton calorimetry" (c.f. iPoh J 1 1 994h PI If proton 
calorimetry is strongly violated, then M82 and NGC 253 and, 
by extension, other starbursts could in fact be weak 7-ray 
sources. 

Although 7-ray emission from M82 and NGC 2 53 has been 
sough t for several years with no suc cess (at GeV. |Cillis et all 
120051 and at TeV, lAharonian etal] 120051: lltoh et all 12007^ 
the launch of Fermi and the advent of powerful VHE 7-ray 
telescopes has led to recen t detections of both starbursts at 
GeV energies (with Fermi; lAbdo et ai1l2010a|) a nd in VHE 
7-rays (M8 2 with VERITAS lAcciari et all 12009; NGC 253 
with HESS. lAcero et al.l 120091) . These GeV and TeV detec- 
tions constrain the cosmic ray (CR) population in these dense 
star-forming environments. 

In this paper, we discuss the implications of the 7-ray 
detections of M82 and NGC 253. The ratio of the 7-ray lumi- 
nosities to the bolometric luminosities informs the question of 
whether or not these systems are proton calorimeters (TQW). 
The emission also has implications for th e energy density of 
CRs in starbursts (e.g.. lAkyuz et al.ll 1 99 lh . Finally, combined 
with the radio emission, the energy lo sses of CR electrons 
and positrons are constrained (c.f., iPaglione et al.l 



Dom ingo-Santamarfa & Torres! 120051; Persic et al J 



1996 



2008 



Ide Cea del Pozo et al.l |2009a; Rep haeli et alJ l2009l) . Pionic 

7 Here, we consider only CR protons with kinetic energy above the thresh- 
old for pion production. 
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7-rays must be accompanied by secondary positrons and 
electrons; the ratio of the power in these expected electrons 
and positrons to the observed radio emission informs us of the 
energy losses of the CR electrons and positrons. In particular, 
we derive in § [5] the expected synchrotron luminosity from 
the pionic luminosity if synchrotron cooling is the dominant 
loss process. 

In {2] we describe the detections of M82 and NGC 253 at 
GeV and TeV energies. We then interpret the detections as 
7-rays from diffuse CR protons in $3] Our interpretation 
includes comparison of the 7-ray luminosities of M82 and 
NGC 253 with their CR luminosities and their IR luminosi- 
ties ( 33.21i, a nd a discussion of the uncertainties in these esti- 
mates ( 33.31 ). We find that a fraction 0.4 and 0.2 of luminosity 
in > GeV CR protons is lost to pion production in M82 and 
NGC 253, respectively. We discuss the implications for our 
estimates mean for pr oton calorimetry in M82 and NGC 253 
are at GeV energies ( 33.41 ). Other possible sources for the ob- 
served 7-rays are considered in $4] The implications of the 
detections of M82 and NGC 253 for the detection of other 
star-forming galaxies, the starburst contribution to the diffuse 
extragalactic 7-ray and neutrino backgrounds, the dynamical 
importance of CRs in starbursts, and for the physics of the 
FIR-radio correlation are described in ^6] We summarize our 
results in $7] 

2. 7-RAY DETECTIONS 

2.1. Fermi and TeV detections 

lAbdo et al] d2010al) reported the detections of M82 (6.8<r) 
and NGC 253 (4.8ct) with the Fermi LAT instrument. At 
energies above a few hundred MeV, the 7-ray spectrum 
of starburst galaxies is expected to be described by a 
power law spectrum with differential photon fluxes N(E) = 
Np(E I Eq)~ t , where r (~ 2) is the photon spectral index (e.g., 
Paglio ne et alJI 1 996b . Using the GeV data point as a normal- 
ization, we adopt GeV differential fluxes of 



GeV 



phcirT 2 s _1 GeV" 1 (1) 



JVngC253 «0.9^ X 10- 9 (JL \ - „„ - r,„v-l 



ph cm -2 s _1 GeV" 1 .(2) 



The lAbdo etail dToiOa) maximum likelihood analysis of the 
sources find T of 2.2 ± 0.2 ± 0.05 for M82 and 1 .95 ± 0.4 ± 
0.05 for NGC 253; however, the Fermi photon statistics are 
insufficient for ac curate measuremen ts of Y. Using the best- 
fit power laws in lAbdo et all d2010al) gives GeV normaliza- 
tions that are 60% - 70% of those in equations Q]|2] although 
these fits are influenced by the nondetections at lower and 
higher energies. Fermi has not detected either starburst above 
~ 20 GeV, and NGC 253 is undetected below ~ 500 MeV. 
The reported spectra are shown in Figure [T] 

Assuming a power law spectrum from GeV to VHE ener- 
gies, we can combine the Fermi detections with the HESS and 
VERITAS measurements to derive the spectral slope over this 
energy range. The integrated fluxes from M82 reported by 
VERITAS (at 1.3 - 3.8 TeV ^ correspond to Ge V- VHE spec- 
tral slopes r of 2.19 - 2.25 (lAcciari et alJl2009l) . in excellent 
agreement with the measured spectral slope from Fermi. The 
HESS detection of NGC 253 at 220 GeV implies Y = 2.3 
dAcero et al.l 120091) . steeper than the best-fit photon index 
from the Fermi detections, but within the quoted errors. We 
adopt T = 2.2 for M82 and T = 2.3 for NGC 253 throughout 



the rest of this paper. 

Our values of Y are only appropriate if the spectrum is truly 
a single power law between GeV and TeV energies. A spec- 
tral bump at GeV energies will cause an underlying pionic 
power-law spectrum to appear steeper than it really is; con- 
versely, a spectral bump at TeV energies will cause it to ap- 
pear flatter. The possibility of a "TeV excess" is notable par- 
ticularly because such an excess is seen in th e M ilky Way 
dProdanovic et alj|2007t lAbdo et alj2008t see §1431 

FigureQ]shows that the 7-ray spectra of both M82 and NGC 
253 are in reasonable agree ment with most of the previous 



detail e d model pre d ictions (Domingo-Santamarfa& Torres 
1 20051; iPersic et alj 120081; Ide Cea del Pozo et alj 12009a; 
Reph aeli et al.l I20091) . However, these models slightly 
overpredict both the GeV flux and the ~ 400 MeV flux 
of NG C 253 by a factor of - 2-3. The iPaglione et all 
(1996) model greatly overestimates the < GeV flux by a 
factor of ~ 5 (although it used a higher CR acceleration 
efficiency), and predicts a much softer spectrum (r ps 2.7) 
than observed. Overall, the general agreement between the 
theory and observations of 7-rays from M82 and NG C 253 
is encouraging (see also Ide Cea del Pozo et al.l l2009bl) . and 
suggests that m odels of other starbursts, particularly Arp 220 
(Torres 2004 b, LTQ), also predict their 7-ray fluxes with 
fidelity (g6j)- 

2.2. M82 and NGC 253 Gamma-Ray Luminosities 

For a dN/dE = No(E /E m i n )~ r 7-ray spectrum from E„^ n to 
£max, the total luminosity at energies greater than E mm is 



L-y(> E mia ) = 2.4 x 10 39 ergs s" 1 N- 9 D 2 3 5 p. 



7> 



(3) 

where AL 9 =M)/(10 _9 phcm- 2 s" 1 GeV _I ),D 3 .5 =£>/3.5 Mpc, 
and 



(T-2)- 



ln 



1- 



2-r 



(IV 2) 



(r = 2) 



(4) 



Typically, /3 7 « 2— 5. 

The GeV flux reported by Fermi for M82 is AL 9 = 1.9^, 
and the VERITAS detections imply a Y ps 2.2 spectrum ex- 
tending to at least E m . dx = 3.8 TeV. This corresponds to a 
> GeV luminosity of 



L M &2(> GeV) ps 1.9^ x 10™D^ 6 ergs s 



VK)r,2 



(5) 



If we interpret the HESS detection of NGC 253 as part of 
a single r pa 2.3 power law extending at leas t to E m!lx = 
220 Ge V, and use the N-g = 0.9^ value from lAbdo et alj 
d2009al) . we find that the GeV luminosity of NGC 253 is 

i N GC253(> GeV) « 5.6ll 5 g x 10 39 D 2 5 ergs s" 1 . (6) 

As E max — > 00, the luminosity increases by only ~ 25%. 
In light of the harder Fermi best-fit spectrum to the GeV 
data (r = 1.95 ±0.4 ±0.05), the VHE emission might be a 
differ ent s pectr al component ( perhaps pulsar wind nebulae; 
see § 14.21 and Mannheim et al.ll2010l) if the pionic emission 
falls off between 20 and 200 GeV. However, a harder spec- 
trum with lower E max still implies a similar 7-ray luminosity 
(~ 7— 11 x 10 39 D 2 5 ergs s" 1 ) with those assumptions, and the 
Fermi data are not accurate enough to conclude there is a dis- 
crepancy. 
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FIG. 1. — The 7-ray spectra of M82 and NGC 253 from Fermi (solid triangles: Abdo et al. (2010 ij); open tria ngles: analysis in Appendix IaI an d VER- 
ITAS and HESS ( filled squares). We show several models: LT Q (unsealed to S FR and using Kennicutt (1998) S g : long-dashe d; scaled to I Sanders et al] 
(2003) L-tir and Kennicutt 1 1998) T, g : solid, grey; scaled to Sanders et al. 12003) Ltir and S e in Table [2] soli d, black). Ide Cea del Pozo et alj )2009al) and 
IDoming o-Santama ria & Torresl T2005) for M82 and NGC 253 respectively (dotted), Persic etal] f2008l) for M82 and lRephaeli et alj 42009T) for NGC 253 (dashed), 
IPaglion e et al. 1 1 99fJ) for NGC 253 (long-dashed), and TQW (p = 2.0, 2.2, 2.4, dash-dotted). We plot E 2 times the differential flux at each energy. At VHE ener- 
gies, the main source of variation in the models is the CR injection slope p. Note that the models of LTQ are proton calorimetric, and TQW explicitly assumes 
proton calorimetry (see Fig. [2). 



3. INTERPRETATION AS PIONIC EMISSION 

3.1. Motivation for Proton Calorimetry 

The hypothesis that starbursts are proton calorimeters is 
motivated by the short pionic energy loss time in their dense 
interstellar media, 

f ^ 2xl0 V(-^)"\ (7) 
V250 cm J / 

where « e ff is the average d ensity of the ISM encountered by 
CR protons before escape (Mannh eim & Schlickeiserll 19941) . 
and the mean gas density in the starbursts of M82 and NGC 
253 is a few hundred cm" 3 . 

If t v is less than the escape timescale, then the system is a 
proton calorimeter. CRs may escape by advection in galactic 
winds or by diffusion. The wind advection time is 

fwind ~ h/v ~ 2 x 10 s yr h m v~ 5 } )07 (8) 

for a scale height of h = 100/iioopc and wind speed of 
v = vsookm s" 1 . We expect the diffusive escape time to 
go as tw(E) = t (E/E )~^ 2 , where t = 26 Myr at E = 
3 GeV in our Galaxy f rom radioactive isotopes in cosmi c 
rays (cosmic ray clocks; IConneil|[l99"8t I Webber et ail 12003b . 
Little is known ab out the diffusion escape time in s tar- 
bursts (for example. IDomingo-Santamarfa & Torr es 2005 and 
Ide Cea del Pozo et al.ll2009al assume ~ 1 - 10 Myr). If we 
assume that CRs stream out of the starbursts at the average 
Alfven speed v A = Bj^JA-np dKulsrud & Pearcdfl969l) . then 

fdiff « 3 .5 Myr ftioo^oo n 25o > ( 9 ) 
where B200 =.6/(200 and «250 = n/(250 cm" 3 ). 



From equations (jTJ - ©, pionic losses dominate advection 
losses if n e ff > 250 cm" 3 /i7 ( ) v 5 oo and diffusive losses at a few 

GeV if n e ff ^ 40 cm^Bjoo^ioo 3 ' a g am assuming CRs stream 
out of the starbursts at the Alfven speed. We see that the ad- 
vective and pionic loss times are roughly equal in M82 and 
NGC 253, implying that inelastic proton-proton collisions are 
an important loss process for GeV protons, and suggesting 
proton calorimetry (LW06, TQW). 

There are two ways to avoid this conclusion: more effi- 
cient escape (shorter f a( j v or fdiff), or less efficient pion losses 
(longer t n , smaller n e ff)- Both are possible. In particular, there 
may be a fa st wind component with v 1000-2000 km s" 1 
from M82 dStrickland & Heckmanl 120091) . In addition, the 
pionic losses are less efficient if the CRs mainly travel 
through low density gas. The ISM in starbursts is clumpy, 
with most of the volume contained in a low density phases 
(e.g.Xord et alJll996HMao et al.ll2000t [Rodrigue z-Rico et alj 
120041 |Westmoquette et al. 2009). Indeed, there is 7-ray and 
radio evidence that CRs do not penetrate deep into molecular 
clouds in the Ga lactic Center, so that « e ff («) and proton 
calorimetry fails dCrocker et al. 2010a.Q). 

If either advective losses or pionic losses dominate, the 
CR proton and hadronic 7-ray spectra should both be rela- 
tively hard with T w 2.0-2.4 for standard injection spectra. 
This is because the equilibrium one-zone CR proton spec- 
trum is roughly N(E) w Q(E)t(E), where Q(E) oc E~ p is the 
proton injection spectrum with 2.0 < p < 2.4 expected, and 
t(E) is the CR proton lifetime including escape and (catas- 
trophic) pion losses. Both the advective and pionic lifetimes 
are roughly energy-independent, so they preserve the hard in- 
jection spectrum. By contrast, in the Milky Way CR proton 
lifetimes are determined by diffusive escape C^diff E~ 1 ' 2 ), 
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so the resulting GeV to Pe V proton spectra go as E~ 2J (e.g., 
iGinzburg & Ptus kin 1976). The harder 7-ray spectra of M82 
and NGC 253 support either proton calorimetry or strong ad- 
vective losses. 

3.2. Measuring The Calorimetry Fraction 

We define L ff as the power of the starburst in all pion end- 
products, including hadronic 7-rays, neutrinos, and secondary 
electrons and positrons. The ratio of to the injected CR 
luminosity (.Lcr) for CR protons with kinetic energy per par- 
ticle (K) above the pion-production threshold CKth) is then the 
"calorimetry fraction," which measures the degree to which 
starbursts are calorimeters: 



-Peal = 



' -77 



Lcr(K > K t h) 



hife 



(10) 



The lifetime fuf e of CR protons with K > includes all 
losses — pionic, ionization, diffusive, and advective. If pro- 
ton calorimetry holds, then fijf e w t^, « Lqr(K > A"thX and 
^cai — > 1 In what follows, we restrict the energy range over 
which we estimate F ca i to be > 1 GeV. In particular, we use 
observed 7-rays with energies > 1 GeV to estimate L ff for 
CRs with K > 1 GeV. While a significant fraction of both the 
pionic 7-rays and the CRs have energies below 1 GeV, lep- 
tonic emission is expected to become increasingly important 
at lower energies, contaminating the estimate of the proton 
calorimetry fraction. 

The total injecte d CR power 'Lcr l ikely scales with the star- 
formation rate (c.f. lAbdo et al.l2010fl) . and concomitantly, the 
supernova rate, Tsn- Assuming that with each supernova, a 
fraction rj' of its kinetic energy goes to primary CR protons 
with A > lGeV, 

L CR (> GeV) = 3.2 x 10 41 ergs s" 1 £51770.1^,0.1 (11) 

where rj' oi = rf /OA, £51 is the energy of the supernova in 
10 51 ergsi and IWi = r SN /0.1 yr" 1 . 

Even in the proton calorimetric limit (£ oa i 1), the GeV 
7-ray luminosity will be significantly smaller than Lcr- First, 
only ~ 1 /3 of L n ends up as 7-rays; « 3L 7 . Second, a frac- 
tion fin of the pionic 7-rays from CR protons with K > GeV 
will have energies > GeV. We calculate j3 n using the GAL- 
PRQF pionic cross sections dMoskalenko & Strong! 119981; 
Strong & Moskalenkol ll998t IStrong et al.ll2000t) based on the 
work of Der mer ( 1986) (see also | Steckerll970l;lBad hwar et al. 
fl^lSteph ens & Badhwarll 19811) and a K p spectrum from 
1 GeV to 1 PeV, and find that it ranges from 0.9 (p = 2.0) to 
0.5 (p = 2.5). Using these factors, (see eq. [3j 



F c . di a 0.023Dl 5 ^V_ 9 /3 7 /3; 1 £5 1 1 ?7 : 1 1 r 



SN.0.1 • 



(12) 



We expect that Tsn is proportional to the luminosity from 
young massive stars and the star formation rate (SFR). Be- 
cause most of the stellar luminosity is converted into infrared 
light by dust in starbursts, if proton calorimetry holds, then 
the 7-ray flux of M82 and NGC 253 should simply be a frac- 
tion of the total FIR flux. Furthermore, if rf is constant for 
all starbursts, then this ratio of observed fluxes will be con- 
stant in the calorimeter limit, so that starbursts should lie on 

8 Note that equation 1101 ignores the energy dependence of the losses, 
aside from the rest riction that K > K t ^ ; ionization losses should be subdom- 
inant for K > K t ^ (Torres 2004b), while both advective and pion losses are 
roughly independent of energy. If diffusive losses dominate in starbursts, they 
may be more effective at higher energy, as in the Galaxy. 

9 GALPROP is available at http://galprop.stanford.edu 



a linear "FIR-7-ray correlation" in analogy with the FIR- 
radio correlation (TQW). Observations of normal and star- 
burst galaxies do suggest some kind of correlation between 
SFR and L 7 , but this correlation increases faster than lin- 
early (L 7 (> 100 MeV) oc SFR 1 4±03 ; lAbdo et al J 120101 . as 
expecte d if escape is more efficient in low luminosity galax- 
ies (c.f. lStronget"al]l20T0h . 

Following TQW, we assume that SFR is related to the total 
FIR luminosity (L TIR [8- 1000] fim) by L pt = eSFRc 2 , wh ere 
e is an IMF-dependent constant (see, e. g.. iKennicuttl 1 99*81) . In 
the calorimeter limit, 



cal 

GeV-TIR 



'^, 3 . 1>lrt „ 4l »„ 



fin Ljir 



(13) 



where *f?n = (rsN/e)/17Mg depends very modestly on the 
star formation history of the galaxy considered; it varies from 
~ 15 to ~ 23 for continuous star for mation over timescales of 

3 x 10 7 - 10 9 yr dLeitherer et alJ[T999l T QW). 

We scale the TIR luminosities from Sanders e t alj d2003l) 
(see Tables Q] & [2]) to the same distances as the 7-ray lumi- 
nosities in equations ((5]) - ©. For M82 (/3 W = 0.7), we find 
that 



^cT 2 = ^-TIR/eGeV-TlR « 0.4 (E 5 , ^ , * ^ ■ 

For NGC 253, we find that F ca i is 



(14) 



7 NGC253 



£Se™/^v-Tm * 0.2(E 5l ^ l7 )- 1 (15) 



for T = 2.3 (to 220 GeV) and j3„ w 0.6. The flux uncertain- 
ties in the Fermi and HESS detections of NGC 253 indicate 
that the uncertainty in £ c ^p C253 can be significantly reduced by 
more Fermi data. Note that we have used the entire TIR flux of 
these galaxies, while the starburst cores (and not the outlying 
disks) probably domin ate the 7-ray emission (see the treat- 
ments of NGC 253 by iDomingo-Santamari a & Torres! 120051 ; 
iRephaeli et aUl2009l) . 

Alternatively, one may estimate F ca i with a distance- 
independent supernova rate (e.g., from radio source counts). 
In this case, F ca i (eq. |T0b retains its strong distance depen- 
dence. For M82 (D = 3.6Mpc, N-g = 2,T = p = 2.2, ^ = 0.7, 
/3 T = 3.96), 



rM82 
^ Ca l : 



:03(E 51 rii 1 r 1 Dt,r- 1 



3.6 1 SN.O.Ij 



(16) 



where we have scaled to a value of Tsn typically quoted in the 
literature. Similarly, 



pNGC 253 
c a l 



; 0.1 (E 5 \r)' Ql y x D\ 5 T 



-1 

SN.0.1 1 



(17) 



for our adopted NGC 253 spectrum (see eq. IT5T>. 

The higher numbers for F ca \ in equations ([Pil l & ( fl5l ) with 
respect to equations ( TTST ) & ( TFTb are easy to understand. For 
= 1, and Ltir gives us Tsn ~ 0.059 yr" 1 and Tsn ~ 
0.049 yr -1 for M82 and NGC 253, respectively. Thus, the 
nominal values for Tsn = 0.1 yr" 1 in equations ( fTSI l & ( TT7b , 
while wel l in the range of supernova rates quoted for both 
systems ( 33.31) . are larger than those inferred from the total 
FIR luminosity by a factor of ~ 1.7-2. 

Note that our estimates of F c - d i in equations $15[ and ( TTTb for 
NGC 253 are still ~ 2-5 times higher than the HESS estimate 
of ~ 0.05, even though we us e similar or higher SN rates. 
The main reason for this is that lAcero et al.l (120091) assume a 
GeV-to-TeV spectral slope of 2.1, whereas we use a GeV-to- 
TeV spectral slope of 2.3, derived from the Fermi data. With 
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this spectral slope lAcero et all ( 120091) effectively underesti- 
mate the GeV 7-ray luminosity of NGC 253 by a factor of 
~ 3 from 1-220 GeV. 

3.3. Primary Uncertainties in F ca \ 

Other j-ray sources - Any 7-ray source besides pionic emis- 
sion from CR protons lowers our estimate for F c . d \. Although 
it is in principle possible that other sources dominate, e.g., 
the TeV emission, it is likely that the GeV emission is in fact 
pionic. See §H] 

Other IR sources - For our main estimates of F zd \ in equa- 
tions ( fT4b and ( fTBT ). we have used the total infrared light of 
each galaxy to measure the star-formation rates of the 7-ray 
emitting starbursts. However, in NGC 253 only about half 
of th e IR emission comes from the starburst (Melo et al. 
2002). It is also possible that cirrus emission from old stars 
contributes to the observed infrared emission, although we 
do not expect this to be significant within the starburst itself. 
Excluding this additional IR light increases the estimates of 

-Fcal- 

Acceleration efficiency - Higher rf l owers our estimated F ca ] . 
In principle, 7/ can be as high as ~ 1 (Ellison & Eichle rll984t 
lEllison et al.ll20041) . Efficiencies 77' > 1 are also possible if 
additional CR power comes from sources other than SNe. 
We have scaled the above estimates for F ca i using rj = 0.1 
based on our work on FIR-radio correlation (LTQ, TQW), 
which constrains E$\rf to be ~ 0. 1, depending on t he CR 
proton injection spectrum (see also lTorres et alJl2003b iTorresI 
2004b). We emphasize that 7/ is the energy per SN explosion 
in CR protons with energies > 1 GeV and does not include 
low-energy CRs. 

Role of supernovae - Although we have assumed in equa- 
tions dl6b and ( fTTI i (and implicitly assumed in our definition 
of rf) that SNe are resp onsib l e for all of the CRs, this 
has not been settled (see iButtl 120 09). There is evidence 
now that SN remnants accelerate some CRs dTavani et alj 
2010), but other sources may also contribute CRs. The 
7-ray detections of M82 and NGC 253, combined with 
the 7-ray detections of quiescent star-forming galaxies, are 
evidence that 7-ray emission scales with star-formation rate 
dAbdo et al. 2010f|). Ho wever, other possible sou rces of CRs 
includ e stellar winds dOuataert & Loebl 12005b. superbub- 
bles dHigdon & Lingenfelteri 12005b ISeaauist & Stankovid 
1 20071: Butt & Bvkovl 120081) p ulsars dArons & Tavani|ll994t 
Bednarek & Protherod 119971: iBednarek & Bartosikl 120041) . 
and gamma-ray bursts dWaxmanll 19951) . which all presumably 
scale with star formation rate. It is also possible that the 
efficiency of some mechanisms, such as superbubble acceler- 
ation, are different in starbursts. 

Supernova rates - Even if SNe are responsible for CR ac- 
celeration, the SN rates in M82 and NGC 253 are highly 
uncertain. Estimates of Tsn come from stellar popu- 
lation fitting dForster Schreiber et ail 12003b. lin e emission 
(Breg man et al.ll2000t lAlonso-Herrero et al.l l2003). FIR emis- 
sion dMattila & Meikldl2001b . co mparison of radio sources 
with models of SN remnants dvan Buren & Greenhousd 
11994b . and direct searches for SNe dMannucci et all 12003b . 
Methods based on the bolometric emission are complicated 
by the star-formation history, potential IMF variations in star- 
bursts, including the high-mass (> 8 M Q ) slope, the shape of 



the IMF below ~ 1 M Q , and the transition mass between stars 
that do and do not produce SNe. Each of these numbers can 
affect ^ in equation ( fT3l ). although we d o not expect large 
variat ions from galaxy to galaxy (TQW; iPersic &"Re phaeli 
|20"K)I) . 

Methods that use direct searches for SNe or their rem- 
nants are more uncertain, and complicated by biases. For 
example, many of the radio sources identified as SN rem- 
nants in M82 and NGC 25 3 may be compact HII regions 
dSeaquist & Stankovid 120071) . Methods based on the expan- 
sion speed of SN remnants may be complicated by differ- 
ent p hysical conditions in starbursts dChevalier & Franssonl 
12001b . At present, only tw o confirmed SNe have been ob- 
serve d in M82 (S N 2004am: ISinger et alJl2004t iMattila et alj 
I2004t SN 2008iz: iBrunthaler et alj|2009bl 12010b . along with 
three radio transients over the past ~ 30 yr that may be radi o 
SNe dKronberg & SrameM[l985t IMuxlow et al][l99l l2QToh . 
In NGC 253, SN 1940e occured 53." (0.9 g 3 .5 kpc) from the 
galax y's center, outside the starburst itself ( Kow al & Sargentl 
11971b . No radio SN was observed in NGC 253 over 17 years 
of observations, but the implied 95% confidence limits on SN 
rate is weak (< 2.4 vr'^ lLenc & T ingav 2006|). 

Overall, Tsn reported in the literature for M82 and 
NGC 253 span an order of magnitude, from 0.03 yr" 1 to 
0.3 yr" 1 . Early e stimates were very high, with Tsn ~ 0.3 yr" 1 
dRieke et al.ll 1980h . More recent estimates have revised Tsn 
downward to ~ 0.1 vr" 1 (fMuxlow et all 1 1 994t iHuang et all 



1994t Ivan Buren & Greenhousd 11994; Bregman et alj 1 200' 



Forster Schreiber et al 



200lT IMattila & Meikld [200 1 ; 



Lenc & Tingavl 120061; Fenech et alj 120081 120101) . Nonethe 



less, Tsn remains uncertain at the factor of ~ 2 - 3 level, and 
SN rates down to 0.02 yr" 1 a r e possible for both s ystems 
dForster Schreiber et al] 120011: IColina & Perez-Oleal 119921: 
Engelbrac ht et al.lll998b . 

Galaxy distances - The estimates of F C!1 \ in equations ( TBI and 
(15[ do not depend on distance, because they depend only on 
the ratio of TIR and 7-ray fluxes. 

Equations ( fToT ) and dT7l > instead assume Tsn, so that 
the distances do matter. Models of the 7-ray emis- 
sion also often assume some distance an d supernova rate 
(e.g.. rDomingo-Santamaria & Torres! [20051 LTQ). The esti- 
mated distances to NGC 253 vary, from less than 2.3 Mpc 
(Davidge & Pritchetl fl990l) to 3.9 Mpc dKarachentsev et alJ 
2003)13 This range amounts to an uncertainty of a fac- 
tor of ~ 4 in the 7-ray luminosity of NGC 253. Simi- 
larly, distances typically quoted for M82 include 3.3 Mpc 
dFreedman & Madore 1981), 3.6±0.3 Mpc from Cepheids in 
M81 dFreedman et all 1994. and 3.9 ± 0.6 Mpc from the red 
giant branch (Sakai & Madore 1999), amounting to a ~ 40% 
uncertainty in the luminosity of M82. 

3.4. Assessing Proton Calorimetry 

The values of F C!i \ derived in Section[X2]imply efficient pro- 
ton losses in NGC 253 and M82 compared to the Milky Way. 
The 7-ray data imply both systems have F c . d \ ~ 0.2-0.5 (com- 

10 Other estimat es include 2.5-2. 7 Mpc fromjde Vaucouleurs 1 1 978|), 2. 6 
Mpc from Puche & Carignan 1 1988), 2.9 ± 0.5 Mpc from lBlechal 11986]). 
3.3 M pc from lMouhcine et alj (2005), and 3.5 ± 0.2 Mpc from lRekola et alj 
120051). A distance of 2 .5 Mpc to NGC 253 is commonly quoted (e.g. 
Mauersberger et al. 1996), and is used in the Domingo-Santamarfa & Torres 
120051) and Rephaeli et al. 1 200 % models of NGC 253; the HESS analysis 
similarly used 2.6 Mpc lAcero et"al]l2009l) . TQW and LTQ used 3.5 Mpc, 
based on the Hubble Law. 
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FIG. 2. — The ratio of the total pionic 7-ray flux at energies > GeV to the total luminosity from star formation for g alax ies with 7-ray detections or upper limits 
(for Arp 220, see Appendix|A)- See Tables2]&[2] The dashed lines are the calorimetric expectation from equation )13t . scaled to rf = 0.05 and using p n = 0.7 
(p = 2.2). The solid line is the predicted ratio for pionic 7-rays in the fiducial model of LTQ based on the Schmidt Law of star formation and the linearity of the 
FIR-radio correlation. The model becomes calorimetric at high surface densities where the curve flattens. The kink at S s = 0.05 g cirT 2 is a result of the scale 
height changing from 1 kpc for normal galaxies to 100 pc for starbursts. Shading indicates the predicted ratios in all successful models of LTQ (darker: p = 2.2 
only; lighter: all p ) for pionic 7-rays. The open star represents the ratio for NGC 253's starburst core, if it has one half the TIR luminosity of t he entire galaxy 
(Melo et al. 2 002ft. The F SF of M31, LMC, SMC, and the Milky Way are based on their SFRs (Williams 2003; Harris & Zaritsky 2004, 20091; I Yin et"aTfc009t) 
and the_Kennicutt 1998 c onversion factor b etween Lfir and SFR. The F 7 (> GeV) are assumed to be pionic, except for the Milky Way, where the pionic 7-ray 
luminosity comes from Strong et al. 1 2010). 



pared with F ca i < 0.1 in the Milk y Way based on grammage 
estimates and modelling; e.g., |G inzburg & Pt uskinl 119761: 
Garcia-Munoz et alT l987; Engelm ann et alJl990tlJones et al] 
200U iDogiel et al.ll2002t IStrong et alj|2010h . This suggests 
that f,r ~ Escape, even though both galaxies exhibit large-scale 
galactic winds, and that these systems represent the transi- 
tion to proton calorimetry. The case for calorimetry in M82 is 
stronger than in NGC 253, although F cd i is uncertain for NGC 
253 by a factor ~ 2. If we only use the TIR flux o f the core of 
NGC 253 (about half the total; iMelo et alj|2002l) . F ca i would 
be about twice as high, or 0.4 - 0.5 - the same as M82 (open 
star in Fig. [2j. Thus, the small F c . d \ for NGC 253 may sim- 
ply be the result of averaging the 7-rays from the calorimetric 
starburst with the non-calorimetric outlying disk. 

Figure [2] shows the ratio F 7 (> GeV)/i*sF (see Tables Q]& 
|2]), the 7-ray flux above 1 GeV to the bolometric flux Fsf pro- 
duced by young stars, as a function of gas sur face density 
for NGC 253 and M82, as well as for the LMC (Port er etalj 



2009) , SMC (lAbdo et all 1201 Odl) the Galaxy (IStrong et alJ 

2010) , and M31 dAbdo et alj|20 10f), together with upper lim- 
its on M33 from lAbdo et all d2010j) and Arp 220 from our 
own analysis of the Fermi data (see Appendix [Ab. (Note that 
the plotted ratio F 7 (> GeV)/fsF does not include /3 W , since 
we wish to plot observable quantities.) The dashed line indi- 
cates the calorimetric expectation from equation ( fT3l , scaled 
to 77' = 0. 1 and ^ = 0.7, for p = 2.2. 

The solid line is the prediction of the fiducial model of LTQ, 
derived by combining constraints from the Schmidt Law of 
star formation and the observed FIR-radio correlation. At low 
gas surface densities, CR protons easily escape, 7-ray emis- 
sion is weak and Cciv-TiR * s small. However, as the gas sur- 
face density increases, galaxies become more proton calori- 
metric and in sufficiently dense starbursts £cev-TiR asymptotes 
(eq.fTJt. For these galaxies we expect a FIR-7-ray correlation 
(TQW). The discontinuity at S g = 0.05 g cm" 2 is due to the 
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transition in scale height from h=\ kpc (for normal galaxies) 
to h = 100 pc (starbursts) in the LTQ models; in reality, the 
transition between normal galaxies and starbursts is smoother. 
The 7-ray luminosity of the LMC is consistent with the pre- 
dictions of LTQ to within a factor of ~ 2. The standard model 
of LTQ is t uned to reproduce the Milky Way 7-ray luminosity 
given in IStrong et all (gOOO) with S g = 0.01 g cm 2 . A more 
recent estimate revises the Mi lky Way pionic 7 -ray luminos- 
ity down by a factor of ~ 3 (IStrong et al .1120101 : S10); on the 
other hand, the Milky Way gas surface densities compiled in 
lYin et al.l (120091) are also ~ 3 times lower (the star-formation 
rate peaks at ~ 6 kpc, where S ? = 0.003 g cm" 2 ), so that the 
Galaxy is still fairly close to the LTQ prediction. The SMC is 
7-ray dim by a factor of ~ 4, suggesting that CR s escape much 
more easily than expected (Abdo et alj|2010db . Similar be- 
havior is indicated fo r the CR electrons by the radio emission 
of irregular galaxies (Murphv et al. 2008). However, M31 is 
surprisingly 7-ray bright with respect to the LTQ prediction. 
Finally, NGC 253 and M82 appear to be somewhat 7-ray faint 
compared to LTQ's fiducial model, consistent with equations 
tO and (Q3). 

A detection of the ULIRG Arp 220 at the level specified in 
Table [2] would improve our understanding of how F ca \ evolves 
with With average densities in its nuclear starbursts ex- 
ceeding 10 4 cm 3 dDownes & Solomonll 19981) . its pionic loss 
time is less than 10 4 yr (eq. |7), difficult to reach with winds 
(eq. |H). According to the models of LTQ, Arp 220 should 
not have a significantly larger F 7 (> GeV)/fsF than M82 and 
NGC 253. If, instead, Arp 220 is much brighter than M82 
and NGC 253, that implies that either rf is much larger in Arp 
220, or escape is more efficient in M82 and NGC 253 than in 
the fiducial model of LTQ. 

4. OTHER SOURCES OF 7-RAY EMISSION 

Our analysis in §[3] assumes that the 7-ray emission is pri- 
marily pionic. Here, we consider other possibilities that could 
reduce F ca i. 

4.1. Diffuse Leptonic Emission 

Primary and secondary CR electrons and positrons (e ± ) 
also contribute to the 7-ray emissivity of M82 and NGC 253 
via bremsstrahlung and Inverse Compton (IC) of predomi- 
nantly dust-reprocessed starlight. Detailed models of M82 
and NGC 253 indicate that these emission pro cesses are sub- 
domi na nt for energies more than 200 MeV (IPaglione et"afl 
1996; Domingo-Santamarfa & Torres 2005; Persic et al. 
20081 Ide Cea del Pozo et al.1 l2009at lRephaeli e t all 120091: 
de Cea del Pozo et al. 2009b). Even without detailed mod- 



elling, these processes almost certainly do not dominate L 7 in 
the GeV-TeV energy range on energetic grounds. First, more 
energy goes to tt° —> 2j production than to the secondary 
e ± . Second, for a typical ratio of total energ y injected in 
prima ry electrons to protons of 1/50 (e.g., IWarren et al.l 
2005), and even if all the electron energy goes to producing 
7-rays, over 90% (F ca i < 3/50) of the protons would have 
to escape the starburst for the proton contribution to the 
7-ray luminosity to be sub-dominant. Assuming that none 
of the protons interact at all, the overall energy budget of 
the observed 7-ray emission implies that the efficiency of 
primary electron acceleration in M82 and NGC 253 would 
have to be ~ 10 times higher than inferred in the Galaxy. 
Finally, the magnetic energy density (B w 200 pG implying 
Ub ~ 1000 eV cm" 3 ) is predicted to be at least as strong 



as r adiation energy density (£/ ra d ~ 200- 1000 eV cm 3 ; 
e.g., IPaglione et al.ll 1 996t iPersic et al.|[2 008) in many models 
(e.g.~ICondon et alj 1199 It iDomingo- Santama rfa & Torres! 
1 20051 TOW IPersic et alJl2008t Ide Cea del Pozo et al.l l2009a: 
lRephaeli et al.l [2009, LTQ). Thus, synchrotron losses will 
increase the energetics requirements even further. 

However, in specific energy ranges leptonic emis- 
sion can dominate the 7-ray emissivity. Bremsstrahlung 
and IC emission probably make up mo st of the 7- 
ray e mission below 100 MeV ( e.g.. IPaglione et alJ 
1996; Domingo-Santamarfa & To rres! | 2005l: IPersic et al] 
20081 Ide Cea del Pozo et al l I2009ah lRephaeli et alJ 120091: 
Ide Cea del Pozo et al.ll2009 b, LTQ). Bremsstrahlung falls off 
steeply with energy, and is unimportant above ~ GeV. The 
IC spectrum is complicated by the shape of the input photon 
spectrum. Given that the photon SED of starbursts is domi- 
nated by the FIR, for a CR e ± injection spectrum Q(E) oc E~ p , 
we expect the IC photon index to be T\c « p/2+\ at ~ TeV 
energies (Rybicki & Lightman 1979)0 Whether pionic or 
advective losses dominate, the pionic emission will have a 
photon index sw p; thus, IC can dominate the VHE 7-ray 
emission, but only if p is substantially greater than 2.0 in the 
calorimetric or advective limitQ In practice, Klein-Nishina 
effects will suppress the IC luminosity beyond a cutoff 
Ekn ~ 18 TeV (A/80 [mi) for target photons of wavelength 
A. For these reasons, IC is unlikely to contribute significantly 
to the TeV emission unless p > 2.5. 

At pr esent, the Fermi detections presented in Abd o et al.l 
(2010a) only use photons with energies greater than 200 
MeV. Improving the limits on 100 MeV photons is essen- 
tial to determining the leptonic contribution, which should 
begin to dominate at lower energies. Direct detection of 
bremsstrahlung and IC from electrons and positrons would 
have strong implications for the synchrotron radio emission, 
and could test the "high-S g conspiracy" postulated by LTQ to 
explain the radio emission of starburst galaxies (see 

4.2. Discrete ~f-Ray Sources 

Because M82 and NGC 253 are unresolved by Fermi, VER- 
ITAS, and HESS, the 7-ray detections include diffuse emis- 
sion from CRs and emission from discrete sourcesQ The 
high energy particles in these sources responsible for 7-rays 
need not contribute to the general CR population. Many such 
sources in the Galaxy are known to be associated with star- 
formation, in cluding pulsars and a large number of unidenti- 
fied sources dAbdo et al.ll2009aT) . and should be expected in 
abundance in starbursts. 

Because relatively little work has been done on the ex- 
pected properties of such sources in starbursts, it is unclear if 
they could dominate the 7-ray emission from M82 and NGC 
253. A number of star-formation phen omena are known to 
be TeV sources (as reviewed b y, e.g., iGreni er 2008; iHorns! 
120081; iHinton & Hofmannl2q09l). including SN remnants. Pul- 
sar Wind Ne bulae (PWNe) dAbdo et al.ll2009bl) . and possibly 
star clusters (lAharonian e t al. 2007). The observed Galactic 

1 1 Electrons at these high energies, far greater than those observed at GHz 
frequencies, are cooled almost entirely by IC and synchrotron. 

12 For example, if p = 2.2, then TV = 2.2 and Tic =2.1; over 1 dex in 7- 
ray energy, this amounts to only a ~ 30% increase in the ratio of IC to pionic 
7-rays. 

13 Emissi on from AGNs is un likely to be important, because no variability 
is observed l Abdo et al. 2010a) and any AGN luminosity in NGC 253 and 
M82 is small compared to the star-formation luminosity (Brunt haler et al.l 
2009). 
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TeV sources tend to have hard spectra (T w 2.2), similar to 
the observed GeV-to-TeV spectra of M82 and NGC 253. 

As an example, PWNe have a total energy budget set by 
the pulsar rotational energy, E m t ~ 2 x 10 50 ergsP^ 2 m , where 
Pom = P/0.01 s is the pulsar spin period at birth, comparable 
to that injected into CR protons by the SN remnants. The typ- 
ical spindown luminosity is E rot rj 6 x 10 39 Pq 01 Bj 2 ergs s" 1 , 
where B\2 = fi/10 12 G is the pulsar magnetic field strength, 
corresponding to a spindown timescale of ~ 10 3 yr. If all of 
this energy went into 7-ray emission in M82 and NGC 253, a 
pulsar birth rate of ~ 0. 1 yr" 1 could easily power the GeV-TeV 
emission. Howeve r, Galactic 7-ray sources like th e Crab (e.g., 
lAlbert et al]|2008ft. Geminga dYiiksel et alj |2009). and HESS 
J 1825- 137 (Uharonian et alj|2006a!) have GeV-TeV luminosi- 
ties several decades lower than this estimate. For example, the 
total 7-ray luminosity of the Crab (L 7 ps 10 35 ergs s" 1 ) implies 
that ~ 1 4 — 1 5 such objects would be needed to contribute 
significantly to the GeV-TeV emission see n from M82 and 
NGC 253 (see also Mannhei m et all 1201 Oh . Given the pul- 
sar birthrate and spindown timescale, this seems unlikely, but 
we cannot rule out the PWNe in starbursts are much more ra- 
diatively efficient in 7-rays than in the Galaxy, for example, 
through stronger IC losses. 

4.3. A TeV Excess? 

The TeV background of the Milky Way shows a 
"TeV excess" above the expected pio nic background 
dProdanovic et al.1 12007t lAbdo et all 12008ft . Whether it is 
caused by unresolved discrete sources or truly diffuse emis- 
sion is not known; nor is it known whether it is hadronic or 
leptonic. The TeV excess varies with Galactic longitude, be- 
ing strongest in th e Galactic Center and the Cygnus regions 
dAbdo et alj|2008l) . The latitude profile of the Galactic TeV 
7-ray emission supports a hadronic explanation for the TeV 
exces s, but leptonic models are not yet excluded (lAbdo et alJ 
2008). 

The TeV excess is visible in the Galaxy because the pio- 
nic spectrum is s teep; effectively, the Te V excess changes T 
from 2.7 to 2.6 dProdanovic et all 12007ft . If M82 and NGC 
253 have hard pionic 7-ray spectra at TeV energies, then the 
7-rays from ambient CRs interacting with their ISM will bury 
any TeV excess. Furthermore, the simplest explanation for a 
hadronic TeV excess is that some regions of the Galaxy are 
denser and more proton calorimetric, and that the TeV ex- 
cess is simply pionic emission from the normal CR protons. 
This effect is observed in molecular clouds located near CR 
acceleration sites i n the Milky Way dAharonian et alj r2006b: 
lAlbert et al.ll2007t IAharonian et alJl2008ftT but if F ca i rj 1, the 
entire starburst is illuminated this way, and the total pionic 
7-ray luminosity cannot be increased further. 

5. THE GEV-GHZ RATIO: A DIAGNOSTIC OF ELECTRON COOLING 
AND THE FIR-RADIO CORRELATION 

The observed 7-rays from M82 and NGC 253 have im- 
portant implications for the physics of the radio emission of 
starburst galaxies (see LTQ and references therein). If pi- 
onic, L 7 necessarily implies production of secondary e*s, 
which produce synchrotron radiat ion and contribut e to the 
GHz emissivity of the starbursts. Rengarajan] d2005ft . TQW, 
and LTQ have all argued that secondary e dominate the 

synchrotron emissio n from starbur sts. Det ailed models o f 

starburst regions bv iPaglione et al.1 (fl996T). iTorresi d2004b). 
iDomin go-Santamarfa & Torres! d2005l) . iPersic et al.l (2008), 



Ide Cea del Pozo et al.1 d2009aft . and lRephaeli etaT] d2009ft find 
that secondary e ± are the majority of ~ GeV e ± , although the 
number of primary electrons is still within a factor of a few of 
the secondaries. 

The 7-ray to radio ratio provides an important constraint on 
the cooling mechanism of GHz-emitting electrons, if the 7- 
rays are pionic. This ratio can be understood through a simple 
argument as follows. Suppose the protons (and secondary e ± ) 
have an E~ 2 spectrum, with equal energy in each log E bin. 
The protons lose energy to pions; roughly 2 times as much 
energy goes into 7-rays as electrons and positrons. Further- 
more, since vc oc £ 1//2 , the synchrotron emission from each 
log bin in e ± energy is spread over 2 log bins in synchrotron 
frequency. Therefore, if the 7-ray emission is dominated by 
diffuse pionic emission, and if the radio emission is domi- 
nated by secondary e ± , 

vF v {E 1 )=AvF u {E e )f m , (18) 

where £ 7 and E e are the energies of 7-rays and e ± re- 
spectively from CR protons of the same energy, and / syn = 
feynAiife) -1 is the fraction of CR e ± power going into syn- 
chrotron. CR protons produce pionic 7-rays with £ 7 rj §.\K p 
and secondary e with E e sa 0.05^ ss E 1 /2. Since elec- 
trons that emit GHz synchrotron radiation have an energy 
of £ G Hz « 560 MeVBjooV/GHz) 1 / 2 , we compare GeV 7- 
rays and GHz radio emission. The GHz flux of M82 and 
the starburst core of NGC 253 are 9 and 3 Jy respectively 
dWilliams & Bower! 120091) . while the ^F„(GeV) fluxes are 
3 x 10" 12 ergs cm" 2 s" 1 and 2 x 10~ 12 ergs cm" 2 s" 1 respec- 
tively. This implies that f~^ n is ~ 8 for M82 and ~ 17 for 
NGC 253, if all of the 7-ray flux is from their starburst cores. 
If we instead consider the total GHz radio emission of NGC 
253, 6 Jy, then f~^ n w 8. This implies strong non-synchrotron 
losses, consistent with bremsstrahlung and ionization cooling, 
and may be a hint that the "high-E;, conspiracy" advocated 
by LTQ as an explanation for the linear FRC in dense star- 
bursts is operating in M82 and NGC 253. Note that if some 
radio emission were from primaries, this would require even 
greater non-synchrotron losses. If B is much higher than we 
suppose, the e accompanying GeV 7-rays emit at higher fre- 
quencies, but since the radio spectrum only g oes as vF v oc u 3 
dKlein et alJll988HWilriams & Bowe"rj|2009l) . even an order of 
magnitude increase in B only changes our conclusions by a 
factor ~ 2. 

M82 and NGC 253 appear to have steeper GeV-to-TeV 
spectral slopes than p = 2.0, but the basic conclusions are un- 
changed with a more careful analysis. Most of the e ± emitting 
at GHz are expected to be secondaries from CR protons, in- 
jected with a E~ p spectrum. The power going into e ± with 
energy greater than K GHz is L e (K e > K GHz ) = F ca iL C R, P (^ > 
^GHz)/3ir,f(^GHz)/c. where f e is the fraction of pionic lumi- 
nosity going to secondary e , and ^ n>e is the fraction of sec- 
ondary e ± power from CR protons with K p > Kghz that is 
in e with K e > Kchz- The energy going into 7-rays with 
energy greater than a GeV is L 7 (> GeV) = F csl \Lcr, p (K > 
GeV)/3 7r (GeV)/ 7 . Finally the radio luminosity is vh v = L e (> 
■^GHz)(^/GHz) l "''/ 2 /3 syn / sec / syn , where /3 syn is a bolometric 
correction factor and / sec is the fraction of e that are pionic 
secondaries. For p = 2.2, we have from the GALPROP cross- 
sections ^(GeV) = 0.7, PkAKghz) = 0.5, and ft yn w 0.1. 
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Since / 7 w 2/ e , 



/ iy \p/2-i 



L 7 (> GeV) » 30^/ sec ^ ( — 



(19) 



This implies that /"^ is - 5 for M82 and - 7 for NGC 253, 
if all of the 7-ray flux is from their starburst cores. It also 
shows that the simpler estimate in equation ( TT~8T > is a useful 
approximation even when p ^2. 

These estimates are consistent with the idea that most of the 
radio emission in M82 and NGC 253's starburst are from sec- 
ondaries undergoing strong non-synchrotron losses. However 
the exact values of these ratios are still fairly uncertain. The 
main uncertainties are the fraction of 7-rays from diffuse pio- 
nic emission, the fraction of radio emission from primary CR 
e^, the fraction of 7-rays and radio from the starburst cores 
as opposed to the outlying disk galaxies, and the uncertainties 
in the 7-ray fluxes. 

Bremmstrahlung and ionization can easily provide these 
extra losses. For e^ radiating synchrotron at v = z/qhzGHz, 
the densities when the bremsstrahlung and ionization 
cooling timescales are comparable to the synchrotron 
cooling timescale (Wis A synch and f io n/r S yn C h < 1) are 

«eff > 67 cm" 3 z /iwB^n and re BfF > 54 cm" 3 ^ghz^oo, 
respectively (LTQ; [Murp hy! 2009). These losses are also 
suggested by the somewhat flattened GHz synchrotron 
radio spectra observed in starbursts (a f=s 0.7), whereas 
pure synchrotron and IC cooling would lead to steep spec- 
tra with a f=s 1 (iThompson et al.l 120061) . Detailed models 
which do include all of these losses regularly do fit the 
radio spectra of M82 and NGC 253 (iPaglione et all I l996| 
Domingo-S antamari a & Torres! 120051: iPersic et al.l 120081: 

as well 
These 



de Cea del Pozo et al. 2009a; Rephaeli et al. 



as the 7-ray spectra (Ide Cea del Pozo et al.l 



2009 



2009b 



models can improve our interpretation of the relevant loss 
mechanisms. Finally, with more 7-ray data from Fermi, the 
lower energy leptonic 7-ray emission may be detected, which 
would directly constrain the importance of loss processes 
such as bremsstrahlung. 

6. IMPLICATIONS 

A number of previous studies have found implications for 
7-ray bright starbursts. We now discuss these implications in 
light of the current detections of M82 and NGC 253. 

6.1. The Detectability of Other Star-Forming Galaxies 

Our inference of F cal « 0.4 and 0.2 in M82 and NGC 
253 implies that a number of local star-forming and star- 
burst galaxies should be visible with next-generat ion TeV 7- 
ray telescopes like CTA (e.g., lKnodlsed"enl2010l) and with 
additional Fermi data (see also IPavl idou & Fields 1200 It 
ICillis et all 120051 TQW). For a numbe r of galaxies chose n 
from the IRAS Bright Galaxy Survey dSanders et alj 120031) . 
Tables [T]and|2]list the distance, total FIR luminosity and flux, 
estimates of the gas surface density S g , and three determina- 
tions of the 7-ray flux: (1) F^(> GeV), the purely calori- 
metric prediction from equation (|131 l. assuming /3 W = 0.7 and 
* 17 = 1, (2) F| TQ (> GeV), the prediction from LTQ given E g 
and scaled to Ftir (solid line, Fig.[2j, and (3) the observed flux 
F° bs (> GeV) in cases where there exists either a detection or 
an upper limit. Our predictions for F^(> GeV) use equation 
( fT3l with nominal values of /3 V = 0.7 (p = 2.2) and rf = 0.1. 
Tables Q]&[2]provides a useful guide to the detectability of all 



local galaxies to the extent that their TIR and 7-ray light is 
dominated by star formation. 

Table Q] provides results for normal star-forming galaxies 
with low E s that are not expected to be calorimetric. Here, 
F 7 al (> GeV) provides an upper limit to the pionic 7-ray flux, 
and F 7 jT Q(> GeV) provides the prediction based on the nom- 
inal estimate of F(> GeV)/fsF for S g (solid line, Fig. |2]i 
and the star-formation rate listed in the table. Since emis- 
sion from old s tars and di fferent UV opacities can affect the 
TIR emission dBelH l2003h . the 7-ray predictions should be 
scaled to known star-formation rates when possible. Fur- 
thermore, leptonic processes may increase the 7-ray emis- 
sion significantly, especially in low density galaxies. Table 
[2] gives numbers for the starbursts in the IRAS BGS. Here, 
F™\> GeV) « Fjf Q (> GeV)H Uncertainties in Table|2]in- 
clude those in the gas surface density S ? of the starbursts and 
the fraction of TIR light associated with the denser starburst 
regions as opposed to the surrounding galaxy in the whole. 

Note that if they are in fact calorimetric, the starburst/ AGN 
systems NGC 4945 and NGC 1068 in Table |2] should be the 
brightest on the sky a fter M 82 and NGC 253 at GeV We 
note that Abdo et al. (2010b) very recently announced the 
7-ray detection of NGC 4945 (1FGL J1305.4-4928), with 
a flux within 50% of the calorimetric pred iction and very 
near the LTQ prediction; Lenainetal. (2010) has also an- 
nounced the detection of NGC 1068. However, it is still un- 
certain how much of the 7-ray flux from these starbursts is 
from star-formation , and how much is from the Seyfert nuclei 
dLenain et al.ll2010l) . 

The values in Tables Q] and [2] should also give the ap- 
proximate neutrino luminosities of near by starbursts, since 
charged pions decay into neutrinos (e.g., Steckerl fl979l) . In- 
deed, models of M8 2 and NGC 253 predict them to be neu- 
trino sources (e.g.. IPersic etal.1 120081: Ide Cea del Pozo et ail 
2009a; iRephaeli et all 120091) . However. NGC 253 and NGC 
4945 as Southern objects are not detectable with IceCub e 
(though they may be detectable with KM3NET; iKatzl 120061) . 
and M82 is at high declin ation where IceCube's sensitivity is 
weakest (lAbbasi et al.ll2009l) . 

6.2. The Diffuse ^-ray & Neutrino Backgrounds 
from Star Formation 

The detections of NGC 253 and M82 at GeV-TeV en- 
ergies together with our determination of F ca i in these 
systems has immediate implications for the diffuse 7-ray 
and neutrino backgrounds from star formation, as dis- 
cussed in LW06 and TQW ( see al so IPavlidou & Fields! 120021 
iBhattacharva & Sreekumari 120091) . Several recent studies 
present calculations of the star-forming galaxy contribution 
to the 7-ray background, using the 7-ray brightness of nearby 
normal galaxies and starbursts, finding contributions of or- 
der ~ 10-50% fe.g..lField set alJl27nft lMaMvaetalJ l27)Tlt 
IStecker & Ventersl2010l) . Here we present an updated version 
of TQW's calculation, based on the IR background and this 
paper's F cah 

The pionic emission from starburst galaxies should be con- 
centrated above 100 MeV, with a power law spectrum above 
~ 1 GeV to PeV energies if unabsorbed. For an accelera- 
tion efficiency of 77' = 0.1 (see LTQ), for CR protons energies 
larger > GeV, the total integrated 7-ray background above 1 
GeV is, following equation (fTJl i and ignoring absorption and 

14 Note that NGC 5128 (Cen A) has not been included in either Table. 
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TABLE 1 

Non-Calorimetric Galaxies: Predicted & Observed Gamma-Ray Fluxes 













Predicted 


Predicted 














Calculated 


Calorimetric 


T TO 
LI 1^ 


Observed 




Name 


D a 


t b 
i-TIR 


SFR C 




F S a '(> GeV) e 


Fj TQ (> GeV) f 


Fj(> GeV) E 






(Mpc) 


log 10 [L Q ] 


(Moyr 1 ) 


/ — 9 — 1\ 

(ergs cm s ) 


/ — 9 - Is 

(ergs cm s ) 


/ —9 —Is 

(ergs cm s ) 


i —9 —Is 

(ergs cm s ) 


(g cm z ) 


LMC 


0.05 


8.83 


0.2 1 


1.47 X 10~ 5 


3.2 X IQ- 9 ' 


6.2 X 10-"' 


5.7 ± 1.4 X l(r"j 


0.002 k 


SMC 


0.06 


7.86 


0.1' 


5.11 X IO" 6 


1.1 X 10-" 


3.1 X IO- 1 " 


1 x 10- llm 


0.003" 


M31 (NGC224) 


0.79 


9.39 


1.0° 


2.95 x 10~ 7 


6.4 x 10- Ho 


6.3 x 10- 13 ° 


2.6 ± 0.6 x IO -12 ? 


0.0011 


NGC 598 (M33) 


0.84 


9.07 


0.5 r 


1.30 x 10~ 7 


2.8 x IO" 11 


5.5 x IO-' 3 


< 1.5 x 10- 12s 


0.002' 


NGC 6946" 


5.32 


10.16 


2.6 V 


1.69 x 10~ 8 


3.7 x IO" 12 


1.3 x IO-' 3 




0.004' 


NGC 5457 (M101) 


6.70 


10.20 


1.7 W 


6.96 x IO" 9 


1.5 x IO" 12 


2.9 x IO-' 4 




0.002' 


NGC 5194 (M51) x 


8.63 


10.42 


3.6 U 


8.89 x IO" 9 


1.9 x IO" 12 


1.0 x IO-' 3 


<8x 10~ lly 


0.006' 


NGC 3031 (M81) 


3.63 


9.47 


0.3 2 


4.19 x 10- 9 


9.1 x 10- 13 


1.3 x IO-' 4 




0.0015 1 


NGC 3521 


6.84 


9.96 


0.9 V 


3.54 x IO" 9 


7.7 x IO" 13 


2.5 x IO-' 4 




0.0035' 


NGC 5055 


7.96 


10.09 


1.3 V 


3.77 x Ur 9 


8.2 x IO" 13 


2.3 x IO-' 4 


<8x l(r Uy 


0.003' 


NGC 7331 


14.71 


10.58 


5 aa 


4.25 x 10"' 


9.2 x IO" 13 


2.2 X IO-' 4 


<8x 10" 1I >' 


0.0025' 



"Distances from IRAS BGS unless otherwise noted. 
b TIR luminosities from IRAS BGS unless otherwise noted. 

c Th e star-formation rate as inferred from the literature. The TIR flux is likely to be inaccurate as a simple SFR indicator at these luminosities (Bell 

12001 . 

We calculated the bolometric star-formation flux from the SFR, using Fsf = 3.8 X 10A 2 SFR/(4ttZ) 2 ), based on the starburst IR to SFR conversion- 
factor in Kennicutt ( 1998). See cavaets in footnote c. 

e Pionic gamma-r ay fl ux predicted in the explicitly calorimetric limit: F 7 (> GeV) = /J-n-fsF x 1-8 x l(X^(Ejir]' os ^n), using /3 n = 0.7 as a fiducial 
value; see equation [131 

f Pionic gamma-ray flux predicted by the fiducial model of LTQ, using F$p. See solid line in Figure [2] Note that leptonic emission (particularly IC) 
may dominate at the lowest S ? and increase the 7-ray fluxes. 
E Measurement of or upper limit on integrated gamma-ray flux of energies > GeV. 
h Gas surface density. Typical uncertainty in this quantity is ~ 0.3 dex. 

IHarris & Za ritskv 1 200 9]) find that a SFR of 0.2 M yr _1 in the LMC (as used in Fig.ff] 

J Calculated using Porter et al. (2009), assuming that T = 2.7 above 1 GeV. Integrating the total emission in Figure 8 of Abdo et al. (2010e) gives 
similar results. 

Calculated using a total gas mass of 6x 10 8 M Q II srae 1 1997) and R 2 5 ~ 4.9 kpc. 

'Harris & Zaritskv (2004) find an average SFR of 0.1 M0 yr~' in the SMC (as used in Fig. [2) over the past few Gyr, with occasional bursts of 
star-formation mo re recently. 

"'Calculated from Abdo et al. (2010d), by integrating the total power above 1 GeV plotted in their Figure 5. 
"Calculated using a total gas mass of 4.5 X 10 s Mq {Israeli 1997) and R25 ~ 3.0 kpc. 



Williams (2003) find an average SFR of 1 M 



yr 



in M31 (as used in Fig.ffJ. 



PCalc ulated from the M ilky Way-scaled GALPROP model in Fig. 2 of Abdo et al. (2010f). See also the upper limit in Blom et al. 11999). 
IFrom Kennicutt 1 1998). The peak gas surface density compiled in Yin et al. 12009) is comparable. 
IGardan et al. (2007) and references therein find star-formation rates of 0.3 — 0.7 Mq yr -1 in M33. 

s Calculated from Abdo et al. (2010f), scaling from the detection of M31 to the upper limit on M33 for the GALPROP spectral template. 
'From lKeimicuta (W% . 

"This system also has a central dense starburst component with S g ss 0.04 g cnT 2 (Kennicutt 1 998), that may be calorimetric, and amounts to ' 
of the total star f orma ti on rate. 
v Scaled from Lerov et al. (2008) to the distance li sted here. 

"Calculated from Kennicutt et al. (2008) using the Kennicutt 1 1998) Ha luminosity to SFR conversion. 

x This system also has a central dense starburst component with S s ft: 0.06 g cm- 2 (Kennicutt 1998), that may be calorimetric, and amounts to ' 

of the total star formation rate. 

REGRET upper limits from lCillis et"al] <2T)0"5l) . 

z Kennicutt et al. (2008) find an Ha luminosity equivalent to 0.5 Mq yr _I , while Davidge 12006) find a star-formation rate from 10 to 25 Myr ago 
(roughly the typical lifetime of GeV CR protons in Milky Way-like galaxies) of 0.1 Mq yr~' . 

"' Thilker et al. (2007) compare star-formation rates derived through several indicators and find them to be 4.4 — 6.3 Mq yr~' . 



10% 



■ 25' 



redshift effects, 



F 7 (> GeV) w 1.0 x 10" 



rburst real 



F,^F™ GeV s' 



-1 -2 

cm sr 



>„' ^-burst peal p- 1 
V0. liO.75 r 0. 5-^20 



GeV s 



-1 -2 - 

cm sr 



(20) 

for r = 2.2, where f 2 o IR i s the tota l diffuse extragalactic TIR 
background in units of 20 nW m" 2 sr" 1 , ./J "I' = /bm-st/0.75 is 
the fraction of the TIR extragalactic background produced by 
starburst galaxies (see TQW), and F§f = F ca i/0.5 is the av- 
erage calorimetric fraction of these starbursts. As with indi- 
vidual starbursts, the observable neutrino background must be 
comparable to the pionic 7-ray background, and have a simi- 
lar spectrum (c.f. LW06). 

This estimate for the 7-ray and neutrino 
backgrounds from CR protons implies that for 
£ 7 > 1 GeV and V = 2.2, i//^(GeV) - 2 x 



lO-^^/GeVr - 2 ^!/^^^^ 

within a factor of ~ 2 of the current observati ons of the 
extragalactic 7-ray background (lAbdo et al.l 1201 Oct see also 
iKeshet et al.1 l2004i) . This indicates that starburst galaxies 
can be a major source of the 7-ray background. More 
detailed modelling of redshift evolution are necessary to get 
the spectral dependence correct; most of the star-formation 
in the Universe occurs at z ~ 1 (e.g., iHopkins & Beacoml 
2006), so redshift effects will be significant. Furthermore, 
the Universe becomes opaque to 7-rays with observed energy 
> 100 GeV at z = 1 (> 50 GeV at z = 2; > 20 GeV at a ll 
reasonable z; e.g., Gil more et all 120091: iFinke et al.l 120101) : 
the 7-ray background above this energy will cascade down 
to lower energies. The corresponding neutrino background 
will be affected by redshift but not by opacity. At energies 
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TABLE 2 

Possible Calorimetric Galaxies: Predicted, & Observed Gamma-Ray Fluxes 



Predicted Predicted 
Calorimetric LTQ Observed 



Name 


D a 
(Mpc) 


£-riR b 
log lo [L ] 


Ftir° 
(ergs crrT 2 s ) 


p^(> GeV) d 
(ergs cm -2 s ) 


I TO 

F*t (> GeV) e 
(ergs cm -2 s ) 


F 7 (> GeV) f 
(ergs cm -2 s ) 


(g cm 2 ) 


IVIOZ. \ IN VJ V, JU Jt f 


3.63 


10 77 


1 42 X 10~ 7 


1 OS v IfT 11 

J.L/O A 1U 


2 X 10 -11 


1.3 x 10~ 11 


24 h 


i> vJ iJJ 


3 50' 


10 54J 


9 09 x 10~ 8 


1 97 x 10 -11 


1.2 x 10~" 


f\ ^ 4- 1 ^ v in -12 

U.J ZH A 1L/ 


15 k 


l > vJ» . '■t 


3 92 


10 48 


f\ 0% y in -8 

\J.^.J A 1U 


1 y 1fl~ 11 

1 . J J A 1U 


8.5 x 10~ 12 


q 7 J_ q n x 10~ 12 1 


19 m 


wr,r 1 nftS (Mil) 


13 70 


1 1 27 


i iw i n~ 8 

J . 1 J A lu 


U . O^t A 1 U 


1 x 10 -12 n 


■X i n v in~ 12 

J . U HI 1 . U A 1 U 


02 n 


Mr T r CMRI'I 


3 60 


10 10 


injjyi fr 8 


fi f\% v ifr 12 

U.UO A L\J 


6.6 x 10~ 13 


< 4 x io _11 P 


011 ,r 


IC 342 


4 60 


10 17 


2.22 x 10~ 8 


iBTy IfH 2 
t.Oi A LU 


7 2 x 10~ 13 




0.021 


IN 1 tO 


16 47 


1 1 07 


1 37 x 10 ^ 


2 97 x 10~ 12 


1 8 x 10~' 2 


< 4 x iQ-Hp 


14 c| 


NGC 3690/IC 694 


47.74 


11.88 


1.06 x 10~ 8 


2.30 x ur 12 


1.6 x 10~ 12 




2.61 


NGC 1808 


12.61 


10.71 


1.02 x 10~ 8 


2.2i x ur 12 


1.3 x 10~ 12 




0.09'' 


NGC 1365 


17.93 


11.00 


9.87 x lO" 5 


2.14 x ur 12 


1.2 x 10~ 12 


<8x 10- n P 


0.081 


NGC 3256 


35.35 


11.56 


9.21 x lO" 9 


2.oo x ur 12 


1.3 x 10~ 12 




0.281 


NGC 4631 


7.73 


10.22 


8.81 x 10~ 9 


1.91 x 10~ 12 


2.8 x 10-' 3 


<4x Kr u P 


0.021 


Arp220 


79.90 


12.21 


8.06 x 10~ 9 


1.75 x 10~ 12 


1.2 x 10~ 12 


<8x 10-"p 


10 s 


NGC 891 


8.57 


10.27 


8.04 x 10"' 


1.74 x 10~ 12 


1.0 x 10~ 12 




0.081 


NGC 3627' 


10.04 


10.38 


7.55 x lO" 5 


1.64 x 10~ 12 


3.6 x 10~ 13 


<4x 10-"p 


0.041 


NGC 7552 


21.44 


11.03 


7.39 x 10"' 


1.60 x 10~ 12 


6.7 x nr 13 




0.051 


NGC 4736 (M94) 


4.83 


9.73 


7.30 x 10" 9 


1.58 x 10~ 12 


3.5 x 10~ 13 




0.041 


NGC 2903 


8.26 


10.19 


7.20 x lO" 9 


1.56 x 10~ 12 


9.0 x 10~ 13 


<8x 10-"p 


0.081 


ESO 173-G015 


32.44 


11.34 


6.59 x 10- 9 


1.43 x 10- 12 


6.0 x 10- 13 




0.05" 


NGC 660 


12.33 


10.49 


6.45 x 10- 9 


1.40 x 10~ 12 


8.1 x 10~ 13 


<8x 10- n P 


0.081 


NGC 1097 


16.80 


10.71 


5.76 x 10- 9 


1.25 x 10- 12 


7.4 x 10- 13 




O.li 


NGC 3628 v 


10.04 


10.25 


5.59 x Kr 9 


1.21 x 10~ 12 


2.3 x 10- 13 


<4x 10-"p 


0.03 w 


NGC 3079 


18.19 


10.73 


5.15 x 10" 9 


1.12 x 10~ 12 


7.8 x nr 13 


<8x 10 _11 P 


3.71 



"Distances from IRAS BGS unless otherwise noted. 
b TIR luminosities from IRAS BGS unless otherwise noted. 
C TIR flux: Ftir = Ltir/(4ttD 2 ). 

d Pionic gamma-ray flux predi cted in the explicitly calorimetric limit: F 7 (> GeV) = /3-n-FriR x 1.8 x W^ i (E}ir]' Q 05 VI 1 17), using f}^ = 0.7 
as a fiducial value; see equation [T3"1 

e Pionic gamma-ray flux predicted by the fiducial model of LTQ. See solid line Figure[2] In these models, leptonic emission is expected 
to be relatively minor (< 10%) when integrated above GeV energies, although it may comprise up to ~ 25% of the differential emission 
at 1 GeV. 

'Observed gamma-ray flux for energies > GeV, or upper limit. 

g Gas surface density. Typical uncertainty in this quantity is ~ 0.3 dex. 

h We take M s = 2.3 X 1O 8 M JWeifi et alJ2"00"!l) and r = 250 pc for the D adopted: E g = M g /irr 2 . 
'Adopted distance different than in IRAS BGS (3.1 Mpc) for consistency with the rest of this paper. 
J TIR luminosity corrected for larger adopted distance. 

k From Kennicutt 1 1998), but scaled to the CO-tb conversion factor advocated by Mauersberger et al. 1 1996). 

'From the 1FGL source catalog, as announced in Abdo et al. 1 2010b). NGC 4945 is a Seyfert galaxy, and the AGN may contribute some 
7-ray flux. 

m Total gas mass within a radius of 12" (~227 pc at D = 3.92 Mpc) is taken as M s ss 1.7 X 1O 8 M (Maue rsberger et all 19961) . 
ISchinnerer etal] (2000) give a gas mass of M g fa 5.7 X 1O 8 M within r w 1.4kpc, implying S g K. 0.02 g cm" 2 for NGC 1068. 
However, the gas mass is not uniformly distributed in this region, but is concentrated in spiral arms. If we instead use Ej 0.1 g cm~ 2 , 

we find F| TQ (> GeV) = 4.0 X 10~ 12 ergs cnT 2 s~' . 
"Derived from the power la w fit to NGC 1068 as found by Lena in et al.l <20TO|) . 
PEGR ET upper limits fr om lCillis et alHIfJOU) . 
IFrom lRenmcunl fl998T) . 

j M83 has a central starburst region with a central surface density of E s « 0.07 g cm -2 and scale radius of ~ 0.6kpc ILundgren et alj 

12004b . 

s From lDown es & Solomo^ 0213). 

'One member of the Leo Triplet (w ith NGC 3623 and th e starburst NGC 3628). 

"Also known as IRAS 13242-5713. Negishi et al. (2001) give diameter of 1.1 armin, corresponding to ~ 9.45 kpc. Using the K ennicutll 
1 1998) relation between FIR luminosity and star formation rate, we derive ~ 37.7 M yr~' and a surface density of star formation of 
f» 0.54 M yr~' kpc -2 . Assuming the Schmidt Law gives an estimate of the gas surface density of E ? S3 0.05 g cm~ 2 . 

v One member of the Leo Triplet (with NGC 3623 and the starburst NGC 3627). 

llsraell <2009T) gives M g a 1 .5 X 10 8 M Q in the inner 0.6 kpc, implying E ? ~ 0.03 g cm" 2 . On larger scales. lTrwin & Sofuei fl996l) derive 
M g 1.7 X 1O 9 M in the inner r ft! 1.95 kpc, implying again that Eg ~ 0.03 g cm -2 . 
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below 1 GeV, the pionic emission should decline because of 
the decreasing pion production cross-section, as has been 
discussed in previous wo rks (e.g. JProdanovic & Fieldsll2004t 
IStecker & Ventersll2010l) . 

The primary uncertainties are still fw st and F cs \ for star- 
bursts at high z- Importantly, Daddi et alj (1201 Oh recently pre- 
sented CO luminosities of near-infrared selected BzK galax- 
ies at z ~ 1.5. They derive gas masses for these relatively 
normal star forming galaxies of ~ 10 11 M Q and radii of 
R ~ 3 - 6 kpc. These numbers imply gas surface densities 
of E g w 0.3Mn/?J 2 g cm" 2 . Comparing this with Figure|2]we 
expect that BzK galaxies are as calorimetric as M82. These 
galaxies are an important contributor to the total star for- 
mation budget of the universe in the critical redshift range 
Z ~ 1-2, thus strengthening the case for a half-calorimetric 
background as in equation (l20l . However, we emphasize that 
this estimate for the BzK galaxies has significant uncertainties 
(e.g., the CO-to-H2 conversion factor). 

6.3. The Dynamical Importance of Cosmic Rays in Starbursts 

CRs are dy namically important w ith respect to gravity 
in the Galaxy (Boulares & Cox|[r990h . They have recently 
been claimed to be sub-dominant with respect to gravity in 
starbursts because of strong pion losses (LTQ). However, 
CRs may be impor tant in driving w inds in such systems 
(ISocrates et al.l2008h. as in the Galaxy (IChevalier & Fr ansson 
U984tlEverett et al.ll2 008). The observed 7-ray emission from 
M82 and NGC 253 can be converted into a constraint on the 
product of n e ff (see eq. |7]i and the energy density of the CRs, 
Ucr, and hence inform the question of whether or not CRs are 
dynamically important in these systems; 

M*cr > GeV) « ObtV/aevAr, (21) 

where L n (K CR > GeV) sa 3/3; 1 L 7 (> GeV), V is the starburst 
volume, and /c e v is the fraction of the CR energy density in 
CRs with energy above 1 GeV. Taking values for the radius 
and scale height of R250 = R /250 pc and h 100 = h/ 1 00 pc, 

t/ C R"eff « 6100 eV cm-VceV^Vloo^L/V^- (22) 

If we assume R250 = /zioo = 1 and that CRs interact with an av- 
erage density of ~ 250 cm" 3 , t hen Ucr is ~ 3 00/n' y eV cm" 3 
for M82 (comparable with lAcciari et al.l (120091) ) and ~ 
100/g'v eV cm" 3 for NGC 25 3 (extrapolating to 1.3 TeV, 
12eV cm" 3 , about twice that of lAcero et al.l 120091: however, 
they quote n e ff sa 600 cm" 3 ). Fo r a K~ 23 CR s pectrum stretch- 
ing from 10 MeV to infinity (iTorresi l2004bT) . / Ge v ~ 0.25, 
mostly in low (<C GeV) energy CRs. 

There is a degeneracy between Ucr and n e g: at fixed L 7 , 
a small n e ff can be accommodated by having a high Uqr, 
and vice versa. The value of « e g is not obvious since the 
ISM of starbursts is highly turbulent and clumpy, with most 
of the volume filled with gas that is underdense with respect 
to the mean density. The pressure required for hydrostatic 
balance is Phydi-o ~ TrGEgEtot, where E tot is the total sur- 
face density and is approximately equal to E g E3 We find 
that the CR pressure Pcr = Uqr/ '3 is dynamically unimpor- 
tant: PcR/Phydm ~ 0-02/q' v . Alternatively, if we assume that 
■Pcr ~ Phydm, we find that n eff « 0.02/^ v (n) w 5 -20 cm" 3 , 

15 Although the thermal pressure within M82 is an order of magnitude 
less than Phydro (Lord et al. 1996; Smith et al. 2006), the turbulent pressure is 
comparable to Phyjm jSm ith et al. 20061). 



implying that t v w 5 Myr, approximately 25 times longer 
than the nominal wind escape timescale (see eq. |8). This re- 
quires far more efficient CR acceleration than 77' sa 0.1 (see 
§ 13.2b . which we consider unlikely. Therefore we conclude 
that P C r <P hy dro- 

7. CONCLUSION 

M82 and NGC 253 have now been detected in GeV and 
TeV 7-rays, with fluxes roughly comparable to previous de- 
tailed predictions. We have shown that the observed 7-ray 
fluxes imply that a fraction F ca \ sa 0.2-0.4 of the energy in- 
jected into high energy CR protons is lost to inelastic colli- 
sions (pion production) with protons in the ISM (for r)' = 0.1). 
However, F ca \ in the range of 0.1 - 1 can be accommodated 
with different SN e rat es and acceleration efficiencies (see the 
uncertainties in § !3.3l and @. We f ind a significantly higher 
F ca i for NGC 253 than lAcero etaf] (120091) because NGC 253 
has more GeV emission than they expected. The uncertainty 
in F ca i will decrease significantly with more observations by 
Fermi, HESS, and VERITAS. 

A future test of proton calorimetry in M82 and NGC 253 
would be a 7 -ray detection of a ULIRG like Arp 220 (c.f. 
Torres 2004b). Arp 220 is more likely to be a proton calorime- 
ter than M82 and NGC 253, with its extremely high average 
gas density. If M82 and NGC 253 are not proton calorimeters 
but Arp 220 is, the ratio of Arp 220's pionic luminosity to its 
stellar luminosity will be greater than M82 and NGC 253 - 
it will be brighter in 7-rays than expected (see Figure [2] Ta- 
bles Q]&|2]|. Unfortunately, Arp 220's flux is expected to be 
challenging to detect with Fermi, although upper limits alone 
may be constraining (as in the Appendix). Stacking searches 
of ULIRGs may also prove useful. 

Pionic 7-ray emissi on implies second ary e^ production 
in these starbursts (c.f. [Rengarajan 2005j); from the GHz to 
GeV ratio, we found evidence of significant non-synchrotron 
losses. This suggests that bremsstrahlung and ionization 
are important energy loss me chanisms for CR electrons and 
positrons (c.f. Murphv 2009). This would support the idea 
presented in Thompson et al.1 (120061) that these losses flat- 
ten the GHz radio spectrum of starbursts (§01. It would 
also support the "high-S g conspiracy" suggested by LTQ to 
explain the linearity of the FIR-radio correlation for star- 
bursts, whereby bremsstrahlung, ionization, and IC losses 
suppress the synchrotron radio emission of CR electrons in 
starbursts, but proton calorimetry leads to secondary electrons 
and positrons that boost the radio emission. 

Whatever the underlying physics of 7-ray production in 
M82 and NGC 253 is, the high fluxes of these starbursts sug- 
gest that other starbursts should also be 7-ray sources. We 
compile our predictions in Tables Q] & [2] Considering that 
much of the star formation i n the universe at high-z is in lumi- 
nous i nfrared galaxies (e.g.JElbaz et alll999tlCharv & Elbazl 
l200lHPerez-Gonzalez et al.ll2005l:lMagnelli et al.ll2009l) . star^ 
bursts might make up a signifi cant fraction (~ 1/2) of the en- 
tire 7-ray background (e.g . . iPavlidou & Fiel ds 2002, TQW, 
Bhatt acharva & Sreekumar 2009; %.2i . If the hadronic in- 
terpretation of the 7-ray flux holds, the neutrino background 
should also be large (LW06). 

Finally, the conclusion that M82 has F ca] sa 0.4 and NGC 
253 has F ca i sa 0.2 implies that the pion cooling timescale is 
nearly equal to the wind escape timescale, ~2x 10 5 yr for 
these systems. This, in turn, suggests that the CR protons on 
average interact with ISM near the mean density. If this is cor- 
rect, then the CR pressure is significantly below the pressure 
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TABLE 3 

Fermi-LAT 7-RAY FLUXES OF ARP 220 



Property 


Value 




100 MeV - GeV 


r 23 b 

JV(100 MeV) d 


(1.4 ±0.5) x 10~ 9 phcnT 2 s" 1 

1.7 ±0.4 

0.74 

(16 ± 8) X lO" 9 ph crrT 2 s" 1 GeV" 1 


GeV - 100 GeV 


$35 = 

r 35 f 

N(l GeV)' 


(0.30 ±0.05) x 10-' ph cm- 2 s~ ] 
... g 

2.2 

(1.2 ±0.7) X 10 -9 phcnT 2 s~ ] GeV" 1 


GeV - 100 GeV; T = 2.2 


$35° 

N(l GeV) 1 


(0.20 ±0. 15) x 10-' phcrrr 2 s _I 
1.7 

(0.24 ±0.17) X 10-' ph crrT 2 s" 1 GeV" 1 



a Best-fit integrated flux from 100 MeV to 1 GeV. 
b Best- fit photon index between 100 MeV and 1 GeV. 
c Test statistic for power law fit between 100 MeV and 1 GeV 
for Arp 220; the square root is roughly the signficance of de- 
tection. 

d Best-fit differential flux at 100 MeV. 

e Integrated flux from 1 GeV to 100 GeV. 

f Best-fit photon index between 1 GeV and 100 GeV. 

g The best-fit model for Arp 220 has a I^j of 5.0, which was 

the maximum allowed by our source model file. This is why 

the differential flux seems large when the integral flux is small 

compared to M82 and NGC 253. This value is almost certainly 

spurious, considering the predicted faintness of Arp 220. (See 

Tabled) 

h Test statistic for power law fit between 1 GeV and 100 GeV 
for Arp 220; the square root is roughly the signficance of de- 
tection. 

1 Best-fit differential flux at 1 GeV. 

needed to support each starburst gravitationally, and CRs are NN06GI68G, the David and Lucile Packard Foundation, and 

not on a vera ge dynamically important deep within the star- the Miller Institute for Basic Research in Science, University 

bursts (§ 16.31 ). of California Berkeley. A. L. was supported in part by NSF 

grant AST-0907890, and by NASA grants NNA09DB30A 

, , T _. , _ T . , and NNX08AL43G. EW's research is partially supported by 

We would like to thank Matthew Kistler, Boaz Katz and es- ISp and AEC Thifj research hag made uge of ^ 

pecially John Beacom for many useful discussions. B.C.L. is NASA /IPAC Extragalactic Database (NED) which is oper- 

supported m part by an Elizabeth Clay Howald Presidential ated fe ^ M p lsion Laboratory> California Institute of 

Fellowship. T. A. T. is supported in part by an Alfred P. Sloan t u 1 j «_ * *u xt *• 1 a j 

1 — 1 1 , , . F . , . _ , ■• 1 xTAnA a Technology, under contract with the National Aeronautics and 

Sv^Tnmr 18 m ^ y it a Ja grant Space Administration. 

#NNX10AD01G. E.Q was supported in part by NASA grant F 

APPENDIX 
FERMI DATA ANALYSIS FOR ARP 220 

We followed the procedure of lAbdo et alj (2009a), using the publicly available Fermi data reduction software. The analysis 
is reviewed in the available online documentation^] We downloaded data from the Fermi LAT data serve £] for METs of 
239557417 to 286813463, a total of 18 months. Data from within 20° of each source were downloaded. We created an exposure 
cube for the entire sky for this time range. We first divided the Fermi energy range into two broad bands: a low energy bin for 
100 MeV < E~, < 1 GeV and a high energy bin for 1 GeV < E 1 < 100 GeV. The source region had a radius of 10°. The selection 
was done with gtselect. We then selected high quality events with gtmktime. 

Using the exposure cube, we created an exposure map around each source using gtexpmap. Finally, we could perform an 
unbinned likelihood analysis with gtlike. We modeled all of the sources listed in the 1FGL Fermi source catalog within 15° 
of Arp 220, the extragalactic background, diffuse Galactic emission, and Arp 220 itself. In each energy band, we fit a power 
law to all of the point sources, including Arp 220; both differential flux and integrated flux were considered for Arp 220. We 
used the P6_V3_DIFFUSE response function. The Galactic background was modeled with the gll_iem_v02.fit background and 

l 6 Located at http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ ,7 Available at http://fermi.gsfc.nasa.gov/ssc/data/. 

; the unbinned likelihood tutorial, which we followed, is specifically at 
http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/likelihood_tutorial.html. 
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the extragalactic background was modeled with isotropic_iem_v02.txt, both of which are the default models. We first found a 
preliminary fit using DRMNFB, and then used the results of that fit to converge to our final fit with MINUIT. 

We did not detect Arp 220, as expected. Our results are summarized in Table [3] For the high energy band, we considered both 
models where T for Arp 220 was allowed to vary, and models where it was forced to 2.2. In the former, T always was forced to its 
maximum value of 5.0, possibly because of a dearth of high energy photons. Both variants give similar results for the integrated 
number of photons above 1 GeV, but in the r^5 model, the normalization of the differential flux at 1 GeV is much higher. 
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